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Abstract 

New results for scalar Casimir self-energies arising from interior modes are presented for the three 
integrable tetrahedral cavities. Since the eigenmodes are all known, the energies can be directly 
evaluated by mode summation, with a point-splitting regulator, which amounts to evaluation of 
the cylinder kernel. The correct Weyl divergences, depending on the volume, surface area, and 
the edges, are obtained, which is strong evidence that the counting of modes is correct. Because 
there is no curvature, the finite part of the quantum energy may be unambiguously extracted. Cu- 
bic, rectangular parallelepipedal, triangular prismatic, and spherical geometries are also revisited. 
Dirichlet and Neumann boundary conditions are considered for all geometries. Systematic behav- 
ior of the energy in terms of geometric invariants for these different cavities is explored. Smooth 
interpolation between short and long prisms is further demonstrated. When scaled by the ratio of 
the volume to the surface area, the energies for the tetrahedra and the prisms of maximal isoareal 
quotient lie very close to a universal curve. 
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I. INTRODUCTION 



The concept of Casimir self-energy remains elusive. Since 1948, the year of H. B. G. 
Casimir's seminal paper pQ, what is now called the Casimir effect has captivated many. 
Casimir discovered an attractive quantum vacuum force between uncharged parallel con- 
ducting plates. Yet, while Casimir later extrapolated from this to predict an attractive 
force for a spherical conducting shell |2], Boyer proved the self-stress in that case to be 
instead repulsive, which was an even more unexpected result Since Boyer's formidable 
calculation, many other configurations were examined: cylinders, boxes, wedges, etc. The 
literature abounds with these results; for a review see Ref. However, since there are other 
well-known cases of cavities where the interior modes are known exactly, it is surprising that 
essentially no attention had been paid to these. For example, recently we presented the 
first results for Casimir self-energies for cylinders of equilateral, hemiequilateral, and right 
isosceles triangular cross sections fSj, even though the spectrum is well-known and appears 
in general textbooks [HlEj- Possibly, the reason for this neglect was that only interior modes 
could be included for any of these cases, unlike the case of a circular cylinder, where both 
interior and exterior modes must be included in order to obtain a finite self-energy. However, 
the extensive attention to rectangular cavities puts the lie to this hypothesis P-IT3] . It seems 
not to have been generally appreciated that finite results can be obtained in all these cases 
because there are no curvature divergences for boxes constructed from plane surfaces. 

In this paper, as in Ref. [5], we put aside the serious objection that these self-energies 
may be impossible to observe, even in principle. For example, the positive self-energy of a 
spherical shell is not the negative of the work required to separate two hemispheres, which 
must be positive. 

In Ref. [5j we obtained exact, closed-form results for the three above-mentioned inte- 
grable triangles, both in a plane, and for cylinders with the corresponding cross section, 
for Dirichlet, Neumann, and perfect conducting boundary conditions. The expected Weyl 
divergences related to the area, perimeter, and the corners of the triangles were obtained, 
going a long way toward verifying the counting of modes, which is the most difficult aspect 
of these calculations. Moreover, we were able to successfully interpolate between the results 
for these triangles by using an efficient numerical evaluation, and showed that the energies 
lie on a smooth curve, which was reasonably well-approximated by the result of a proximity 
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force calculation. In this paper we show that the same techniques can be applied to tetra- 
hedral boxes; again, there are exactly three integrable cases, where an explicit spectrum can 
be written down. Again, it is surprising that the Casimir energies for these cases are not 
well-known. The only treatment of a pyramidal box found in the Casimir energy literature 
appears in a relatively unknown work of Ahmedov and Duru [H] , which, however, seems to 
contain a counting error. 

In this paper we present Casimir energy calculations for the three integrable tetrahe- 
dra. For each cavity we consider a massless scalar field subject to Dirichlet and Neumann 
boundary conditions on the surfaces. We regulate the mode summation by temporal point- 
splitting, which amounts to evaluation of what is called the cylinder kernel [15], and extract 
both divergent (as the regulator goes to zero) and finite contributions to the energy. We 
also revisit cubic, rectangular parallelepipedal, triangular prismatic, and spherical geome- 
tries with the same boundary conditions. In the end, we explore the functional behavior of 
the Casimir energies with respect to an appropriately chosen ratio of the cavities' volumes 
and surface areas. We also examine limits as the prism length tends to zero and infinity, 
which correspond to the Casimir parallel plate and infinite cylinder limits, respectively. 

A. Point-splitting regularization 

We regularize our results by temporal point-splitting. As explained in Ref. 0, after a 
Euclidean rotation, we obtain 



where the sum is over the quantum numbers that characterize the eigenvalues, and r is the 
Euclidean time-splitting parameter, supposed to tend to zero at the end of the calculation. 
One recognizes the sum as the traced cylinder kernel [12]. Next, we proceed to re-express 
the sum with Poisson's summation formula. 
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B. Poisson resummation 



Poisson's summation formula allows one to recast a slowly convergent sum into a more 
rapidly convergent sum of its Fourier transform, 



By point-splitting and resumming, we are able to isolate the finite parts, which are Casimir 
self-energies, and the corresponding divergent parts, which are the Weyl terms (see Appendix 
A for more detail). 

II. CASIMIR ENERGIES OF TETRAHEDRA 

The three integrable tetrahedra mentioned above are not recent discoveries. They have, 
in fact, been the subject of a few articles [TBHIS]- However, there appears to be only one 
Casimir energy article concerning one of these tetrahedra, which we denote as the "small" 
tetrahedron [14]. These tetrahedra are integrable in the sense that their eigenvalue spectra 
are known explicitly, and there are no other such tetrahedra. We will successively look at 
the "large," "medium," and "small" tetrahedra, as defined below, and obtain interior scalar 
Casimir energies for Dirichlet and Neumann boundary conditions. Although the exterior 
problems cannot be solved in these cases, the finite parts of the interior energies are well 
defined because the curvature is zero, and hence the second heat kernel coefficient vanishes. 

A. Large Tetrahedron 

The first tetrahedron, sketched in Fig. [ij which we denote "large," is comparatively the 
largest or rather the most symmetrical. One can obtain a medium tetrahedron by bisecting 
a large tetrahedron and idem for the small and medium tetrahedra. One should note that 
the terms "large," "medium," and "small" are merely labels, since one can always rescale 
each tetrahedron independently of the others. The spectrum and complete eigenfunction set 
for the large tetrahedron, as well as those of the other tetrahedra, are known and appear in 
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FIG. 1: Large tetrahedron: —x < z < x and x < y < 2a — x. 



With Dirichlet and Neumann boundary conditions, different constraints are imposed on the 
spectrum, that is, on the ranges of the integers k, m, and n. 

1. Dirichlet BC 

The complete set of modes for Dirichlet boundary conditions is given by the restrictions 
< k < m < n. After extending the sums to all of (fc, m, n)-space, we remove all the 
unphysical cases which are k = 0, m = 0, n = 0, k = m, k = n, and m = n while keeping 
track of the 24 degeneracies and compensating for oversubtractions. Finally, after suitable 
redefinitions of some individual terms, the Dirichlet Casimir energy for the large tetrahedron 
can be defined in terms of the function 



(2.2) 



and written as 




oo 



g{p, q, r) + g{p + 1/2, q + 1/2, r + 1/2) - 6g{p, q, q) (2.3) 



p,q,r=— 



oo 



6 g{p + 1/2, q + 1/2, q + 1/2) + 8 g{V3p/2, 0, 0) 



+ 3g{p,0,0) , 
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where the sums extend over all positive and negative integers including zero. (In the third 
and fourth terms only p and q are summed over, while in the last two terms only p is 
summed.) Note that the time-splitting has automatically regularized the sums, and it is easy 
to extract the finite part (the Epstein zeta functions Z^, Z^h, etc. are defined in Appendix 
B), 

^^""^ = -1-^ [^3(2; 1, 1, 1) + Zsbi2; 1, 1, 1)] + -^C(3/2)L_8(3/2) + -Lz,,(3/2; 2, 1) 



96 72 

where (the prime means the origin is excluded) 

oo 

5^' {m' + 2nY' = 2as)L^sis). (2.5) 

m,n=— oo 

The energy then evaluates numerically to 

EiD)_ 0-0468804266 ^ ^^.6) 

The divergent parts, also extracted from the regularization procedure, follow the expected 
form of Weyl's law with the quartic divergence associated with the volume V, the cubic 
divergence associated with the surface area S, and the quadratic divergence matched with 
the edge coefficient 

(fl) ^ 3V _ C 
Here and subsequently, the edge coefficient C for a polyhedron is defined as [20] 

where the aj are dihedral angles and the Lj are the corresponding edge lengths. The above 
expression for the divergences will be the same for all subsequent cavities with Dirichlet 
boundary conditions. 

2. Neumann BC 

In the case of Neumann boundary conditions, the complete set of mode numbers must 
satisfy < k < m < n, excluding the case when all mode numbers are zero. The Neumann 
Casimir energy can be defined in terms of the preceding Dirichlet result as 



2 C(3/2)L_8(3/2) + ^25(3/2; 2, 1)J , (2.9) 
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which gives us a numerical value of 



E 



(N) 



0.1964621484 



(2.10) 



The divergent parts also match the expected Weyl terms for Neumann boundary conditions. 
We note that the cubic divergence's coefficient changes sign when comparing Dirichlet and 
Neumann divergent parts: 



3V S C 

+ 7^ — ^ + 



(2.11) 



This form is also obtained for all the following calculations involving Neumann boundary 
conditions. 

B. Medium Tetrahedron 




FIG. 2: Medium tetrahedron: < z < x and x < y < 2a — x. 



The eigenvalue spectrum of the medium tetrahedron, shown in Fig. |2| obtained by bi- 
secting the large tetrahedron in the z = plane, is of the same form as that of the large 



tetrahedron [Eq. (2.1)] with different constraints. 



1. Dirichlet BC 



The complete set of mode numbers for the Dirichlet case satisfies the constraints < 
m<n<k<m + n. Following the same regularization procedure used in the preceding 
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cases, we obtain the Dirichlet Casimir energy in terms of the Dirichlet result for the large 
tetrahedron, 



3C(3/2)/3(3/2)-(l + v^)7r2l, (2.12) 



2 ^ 967ra 
where we used [T9] 

oo 

E' im' + = KisMs) . (2.13) 



m,n=— oo 

The Casimir energy evaluates to 

0.0799803933 ^ 

Here the function /3 is also defined in Appendix B. 

2. Neumann BC 

With Neumann boundary conditions, the complete set of mode numbers is restricted to 
0<m<n<k<m + n, excluding the all-null case. As with the Dirichlet case, the 
Neumann Casimir energy for the medium tetrahedron can be expressed in terms of the 
Neumann result for the large tetrahedron: 



1 „(N\ 1 



0.1997008024 



a 



3C(3/2);g(3/2) + (l + v^)7r^ 



(2.15) 



C. Small Tetrahedron 

The small tetrahedron (Fig. |3| may be visualized as the result of a bisection of a medium 
tetrahedron along the plane y = a. The form of the eigenvalue spectrum for the small 
tetrahedron is different from the previous two tetrahedra but the same as the cube's^: 

,2 



^lmn=^{k' + m' + n') . (2.16) 
The Dirichlet case is the aforementioned "pyramidal cavity" considered in Ref. 



^ This spectrum is actually the same as that for the other tetrahedra, given in Eq. (2.1 1, with the additional 
restriction that m + n + k he even. 
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FIG. 3: Small tetrahedron: < z < x and x < y < a. 



1. Dirichlet BC 



The modal restriction for the complete set is < /c < m < n. The finite part obtained is 
thus 



E 



1 



^ ^3(2; 1, 1, 1) + -^C(3/2)L_8(3/2) + ^C(3/2)/3(3/2) 



327T 



64 72 96 



which evaluates to 



E 



0.10054146218 



(2.17) 



(2.18) 



This result differs from that of Ref. [M]. The discrepancy appears to stem from a mode- 
counting error in Ref. jH], and the result found there is likely wrong. 



2. Neumann BC 



For the Neumann case, we again find the same condition that the mode numbers must 
satisfy < k < m < n excluding the origin. The Neumann Casimir energy is derived to be 

1 



with a numerical value of 



E 



(N) 



C(3/2) 2L_8(3/2)+/3(3/2) 



0.2587920021 



(2.19) 



(2.20) 
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III. CASIMIR ENERGIES OF RECTANGULAR PARALLELEPIPEDS 



Amongst the geometries considered for Casimir energy calculations, rectangular paral- 
lelepipeds are the most straightforward. Their eigenfunctions and eigenvalues are well known 
and as such have been subject of many articles p4T3]. We rederive a few of these results 
in the following paragraphs. We consider a generic rectangular parallelepiped of length a, 
height b, and width c. The eigenfunctions for a Dirichlet rectangular parallelepiped are the 
well-known products of three sine functions. The spectrum is the familiar expression 



kmn 



+ 



62 



+ 



n 



(3.1) 



A. Dirichlet BC 



The complete set of mode numbers in the Dirichlet case satisfy the restrictions k > 0, 
m > 0, and n > . The Casimir energy may be written in terms of Epstein zeta functions 
and the ratios: x = iP/'^Y cind a = (c/a)^ 



E 



327r2 



^3(2; 1,X,^) 



X ^2(3/2; l,x) 



647r 
96 V 



/^Z2(3/2;x,a) + v^Z2(3/2;l,a 



1 1 

+ 



(3.2) 



B. Neumann BC 



For Neumann boundary conditions the complete set is given by > 0, m > 0, and n > 0, 
excluding the case where they are all null. In terms of the Dirichlet result we obtain 

1 



E 



(TV) 



E 



327ra 



X^Z2(3/2; X, a) + 7^^2(3/2; 1, a) + v^Z2(3/2; 1, x) 



(3.3) 



IV. CASIMIR ENERGIES OF A CUBE 



The cube is a special parallelepiped with equal length, width, and height. Our results for 
the generic parallelepiped therefore apply for the particular case of a = 6 = c or x = (t = 1. 
This particular geometry has also been the subject of prior inquiries, for example Ref. [8], 
so we are simply rederiving these results. 
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A. Dirichlet BC 



The Dirichlet Casimir energy for a cube of edge length a is simply 

TT 

(4.1) 



-^Cube ~ ^ 



^ Z3(2;l,l,l) + -|-C(3/2)/3(3/2) 



327r2 ■ / 32 

from which we obtain the finite part, a result which matches that of Ref. [TT] : 



(^) _ 0.0157321825 

-^Cube ~ • l^-^J 

Ctr 



B. Neumann BC 

With the same modal restrictions as for the parallelepipedal Neumann cases, the Neu- 
mann result can be related to the Dirichlet result with 

4uL = 4uL-^C(3/2)/3(3/2), (4.3) 
ovra 

which gives a numerical value already confirmed in Ref. [TT] : 

^„_a™4722^ (4.4) 



V. CASIMIR ENERGIES OF TRIANGULAR PRISMS 

Since infinite triangular prisms are soluble cases [5], one would also expect finite prisms 
to be soluble [Ml |2T] . Indeed, one can also find the interior Casimir energies of finite trian- 
gular prisms of right isosceles, equilateral, and hemiequilateral cross-sections. The spectra 
differ slightly from the infinite cases with the replacement of an integral over longitudinal 
wavenumbers by a sum over discrete longitudinal eigenvalues. 



A. Right Isosceles Triangular Prism 



1. Dirichlet BC 



The spectrum for a prism of right isosceles cross-section, illustrated in Fig. |4| is 



>^lmn = ^i^^ + ^^) + -^k\ (5.1) 
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FIG. 4: Right isosceles prism. \DE\ = \EF\ = a, \DF\ = aV2, and \BE\ = b. 



The complete set of eigenmodes for this Dirichlet case is characterized by < n < m and 
< k. The Casimir energy, in terms of x = {b/aY, is, therefore, 

1 



E- 



(D) 
RIsoP 



647^^3(2; 1. 1. X) + ^C(3/2)/;(3/2) + ^Z,(3/2; 1. x) + ^ 



TT 



TT 



2. Neumann BC 



:i + V2) 

96 



Z2(3/2;l,2x) 
(5.2) 



For the Neumann case, the constraint on the mode numbers is again less strict, with 
0<n<m,0<k excluding k = m = n = 0. In terms of the Dirichlet result, we find: 

1 



^RIsoP - ^RIsoP 32^^ 



2C(3/2)/3(3/2) + 7x^2(3/2; 1, x) + 2^Z2{3/2; 1, 2x) 



(5.3) 



B. Equilateral Triangular Prism 



The spectrum for an equilateral prism of height b, shown in Fig. [5} is 

-— ^(m -mn + n) + —k . 
9a^ 0^ 



X 



kmn 



(5.4) 
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FIG. 5: Equilateral prism. \DE\ = \EF\ = \FD\ = a, and \BE\ = b. 



1. Dirichlet BC 



The constraint on m, and n for the complete set of modes is the same as in the Dirichlet 
right isosceles case. The Casimir energy, in terms of x = {b/a)"^, is thus derived as 



E- 



(D) _ 

~ a] 7r2 



1 



[Zs{2; 3, 9, 16x) + Z,,{2; 3, 9, 16x)] + — C(3/2)L_3(3/2) + —^26(3/2; 1, 



247r 



+ ^^^2 3/2; 9, 16y 



(5.5) 



where we used 



{m' + 3n')-' = 2(1 + 2'-'')as)L^3{s) . 



(5.6) 



m,n=—oo 



This particular case was also considered earlier by Ahmedov and Duru [22], although their 
result appears misleading. 



2. Neumann BC 



The Neumann constraint is also the same as that for the Neumann right isosceles tri- 
angular prism. Similarly to previous cases, we relate the Neumann result to the Dirichlet 
result, 

1 



E- 



{N) _ j^{D) 
EqP ~ -^EqP 



24na 



5 C(3/2)L_3(3/2) + 2 ^2^(3/2; 1, 3) + 72^ ^2(3/2; 9, 16x) • (5.7) 
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C. Hemiequilateral Triangular Prism 
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FIG. 6: Hemiequilateral prism. \DF\ = \EF\ = a/2, \DE\ = a, and \BE\ = b. 



The hemiequilateral triangular prism (Fig. |6]) or prism with cross-section being a triangle 
with angles (vr/2, tt/3, tt/6) shares the same spectral form as the equilateral triangular prism, 



x2 _ 167r%_2 

^kmn 



TT 



(m — mn + n ) + —k . 



9a2 ' ' 62 

They differ, however, in the constraints for each boundary condition. 



(5.8) 



1. Dirichlet BC 

The complete set of modes must satisfy < A;, and < ri < m < 2n. Again, in terms of 
X = {b/aY, we find that the Dirichlet Casimir energy is of the form 

El^L - Ie'H + ^Z.(3/2; 3. ^ - ^ (v^ + 4) ■ (6.9) 

2. Neumann BC 

The completeness constraint for the Neumann case is again less strict than for the Dirich- 
let case with 0<k,0<n<'m<2n excluding the origin. In relation to the previous result, 
we write 

<ip - i« - ^z.m, 3. m -^S^s,^). (5,10) 
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VI. CASIMIR ENERGIES OF A SPHERE 



The sphere is also one of the geometries most often the subject of Casimir energy calcula- 
tions (For more complete references see Ref. ^.) We report results found in the literature for 
a sphere of radius a satisfying Dirichlet and Neumann boundary conditions. A noteworthy 
difference between the calculations of the energies for tetrahedra and prisms as compared to 
a sphere is that in the polyhedral cases only the interior modes are considered (the exterior 
modes are unknown) whereas in the spherical case both interior and exterior are (necessarily) 
included to cancel the curvature divergences. 

A. Dirichlet BC 

The Dirichlet Casimir energy of a sphere is well known and may be found in Ref. 

.(D) 
^Sphere 

B. Neumann BC 



The Neumann result is also well known and can be found in Ref. 



_ 0.0028168 

-^Sohere ~ ^ • \^-^) 



4Ze = -^-'^^- (6-2) 



0.223777 



VII. SYSTEMATICS OF CASIMIR ENERGIES 



As indicated in the introduction, the relation between the self-energy of a system and 
its geometry is not obvious. (We set aside the more serious physical difficulty as to the 
meaning and the observability of Casimir self-energies.) Having additional data, such as 
the self-energies of tetrahedra and finite prisms, may help in shedding some light on this 
problem. Our analysis is similar to the one applied earlier to infinite prisms [5]. In terms of 
the volume V and the surface area S of the bodies, the dimensionless scaled Casimir energies, 
Esc = E X V/S, are tabulated in Table [l| and are plotted against the corresponding isoareal 
quotients, Q = SQirV'^/S'^ in Fig. [7| It is also possible to look at the prisms in a more 
revealing light by plotting their scaled energies with respect to the parameters Q and b'^/A 
(Fig. |8|. The corresponding results for Neumann energies are given in Fig. [9| 
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Cavity 


Q 


4? 


4f 


Small T. 


0.22327 


-0.00694 


-0.01787 


Medium T. 


0.22395 


-0.00696 


-0.01739 


Large T. 


0.27768 


-0.00552 


-0.02315 


Cube 


0.52359 


-0.00262 


-0.04755 


Spherical Shell 


1 


0.00093 


-0.07459 



TABLE I: Scaled energies and isoareal quotients. 
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FIG. 7: Scaled Dirichlet energies vs. isoareal quotients. From left to right, the circular markers 
indicate results for the small and medium tetrahedra (which cannot be resolved on this graph), the 
large tetrahedron, the cube, and the sphere. The square, filled triangle, star, and empty triangle 
markers correspond respectively to square prisms (parallelepipeds with u = 1), and equilateral, 
right isosceles, and hemiequilateral triangular prisms, respectively. The prism energies become 
more negative for 6—7-0, less negative for h — )• oo. Note that the cusps, corresponding to the 
maximal isoareal coefficient for a given class of cylinders, also lie close to a universal curve that 
passes through the tetrahedral points. 
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FIG. 8: Scaled Dirichlet energies of prisms vs. isoareal quotients Q and b'^/A, where A is the cross- 
sectional area. Starting from the lowest Q-values, the curves correspond respectively to hemiequi- 
lateral, right isosceles, and equilateral triangular prisms, and square prisms (parallelepipeds with 
a = 1). The square prisms' curve goes through the cube's data point (displayed prominently). 

VIII. LIMITING CASES OF PRISMS 

For the prisms we considered, with square and equilateral, hemiequilateral, and right- 
isosceles triangular cross sections, we can examine two limits which reduce to known cases. 
With b as the height of the prism, the limit 6 — > must coincide with the classic Casimir 
case of parallel plates of area A, 

-A. (8.1) 

And with 6 — > oo we must recover the energy for infinite cylinders given in Ref. [5] , 

E ^ ^cyl6, (8.2) 

in terms of the energy /length for the cylinder, £cy\- For a prism of cross section A and 

cross-sectional perimeter P, the scaled Casimir energy is 

V EAb . , 

= ^ 5 = (2A + P6)' ^^-^^ 

where 5* denotes the total surface area of the prism and V its volume. Thus, in the limit of 

vanishing height, 

b^O: Esc^ TTT. (8.4) 

17 



0.0 



0.2 



0.4 



0.6 



0.8 



1.0 



Q 



FIG. 9: Scaled Neumann energies vs. isoareal quotients. From left to right, the circular markers 
indicate results for the small and medium tetrahedra (which cannot be resolved on this graph), the 
large tetrahedron, the cube, and the sphere. The square, filled triangle, star, and empty triangle 
markers correspond respectively to square prisms (parallelepipeds with a = 1), and equilateral, 
right isosceles, and hemiequilateral triangular prisms, respectively. The energies in this case are 
always negative. Note that the cusps, corresponding to the maximal isoareal coefficient for a given 
class of cylinders, again lie close to a universal curve that passes through the tetrahedral points. 

In the limit of infinite length, 



b^oo: Esc ^ (Scyib) — . 




These limits are exact. 



It may be interesting to express these limits in terms of the isoareal quotients: 



Q 



367ry2 
53 



{2A + hP) 



367r6M2 



3' 




whence, 



TT' 



.3 



6^0: ^ 



640 Q 
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In this limit Q — t- 0. A similar expression may be obtained in the large b limit if we use the 
proximity force approximation as discussed in Ref. [3]: 



ScylA 



TT 



p2 

A 



368640 

which becomes exact as the smallest angle of the triangle vanishes. In this approximation 



6 oo : E; 



Sc 



vr^ 1 



10240 Q' 



(8.9) 



where again Q — )■ 0. In Fig. 10 we show how the long and short prism limit are approached 
by our data. 

Let us examine the 6 — )• limit for the example of the square prism, where the Dirichlet 



energy is given by Eq. (3.2), with a = 1. Using the Euler-Maclaurin summation formula, 



we find the asymptotic limit of that expression for short prisms to be 



E, 



D 



' 90 6^,^, 15 6^ 
1---C 3 + 



b/a 0. 



^.10) 



There are only exponentially small corrections to this result. Similarly in the long distance 



(infinitely long cylinder) limit, b ^ a = c, 



167ra a 



TT 



G-C(3) 



C(3/2)/3(3/2) 



TT 



a/b 0, 



(8.11) 



again, up to exponentially small terms. (Here G is Catalan's constant.) Note that the leading 



terms in the expressions are the correct limiting forms: that in Eq. (8.10) is Casimir's result 
for plates of area a^, and that in Eq. (8.11) is that for a square cylinder found in Eq. (5.4) of 



Ref. [5j. In Fig. 11 we plot, parametrically, these asymptotic limits against the exact form of 
the isoareal coefficient, e.g. Q = 97r6^a/2(a + 26)^ for the square prism. We note that both 
limiting forms are well reproduced by the data for the prisms, all the way down to the cusp 
which occurs for the maximal value of the isoareal coefficient. In fact, the accuracy of the 



asymptotic formulas is remarkable: At the cusp, b = a, Eq. (8.10) gives for the scaled energy 



of the cube the value -0.00261078, while Eq. ( |8.11[ ) gives -0.00261479, both differing by 
less than 0.5% from the exact value -0.00262203. 



IX. CONCLUSIONS 



In this paper, we have extended the work of Ref. [5] from infinite cylinders to finite 
prisms and to integrable tetrahedra. The previous work was essentially two dimensional, so 
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FIG. 10: EQV/S plotted versus Q in the 6/a — )■ (left panel) and 6/0—7-00 (right panel) limits for 
Dirichlet prisms. The curves consisting of square, filled triangle, star, and empty triangle markers 
correspond respectively to prisms of square, and equilateral, right isosceles, and hemiequilateral 
triangular cross-sections. The curves converge to the expected value of -tt^/6A0 = -0.0484473 
as 6/a — )• 0. In the 6/a — )• 00 limit, they converge respectively to the values of 0.00213, 0.00269, 



0.00274, and 0.00277 which are simply obtained from Eq. (8.5). These limits are not particularly 
close to the proximity force approximation value of 7r^/10240 = 0.00303. 

it was possible to give closed form results for the Casimir self energy. This is apparently 
not possible, at least not currently, for the cases considered in this paper. Nevertheless, our 
answers are expressed in terms of Epstein zeta functions and other well-known functions, so 
numerical results of arbitrary accuracy are available. 

The emerging systematics are very intriguing: Not only do the three integrable tetrahedra 
and the cube and the sphere seem to lie very close to a universal curve (there is a very slight 
discrepancy in the case of the small/medium tetrahedra) but the maximal isoareal limits of 
the triangular prisms line up as well. The square prisms are well described by the asymptotic 



formulae (8.10) and (8.11), and similar formulae exist for the other prisms. These results are 



not yet conclusive, since the cases we can evaluate are limited. Numerical work will have to 
be done to explore the geometrical dependence of the Casimir energy of cavities composed 
of fiat surfaces of arbitrary shape. 

Work on other boundary conditions, in particular electromagnetic boundary conditions, 
is currently under way. Unlike for cylinders, the electromagnetic energy of a tetrahedron is 
not merely the sum of Dirichlet and Neumann parts; there is no break-up into TE and TM 
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FIG. 11: Asymptotic fits for square prisms. The short height (6/a — )• 0) (left panel) and long 



height (6/a — )• oo) (right panel) asymptotic energies are shown, corresponding to Eqs. (8.10) and 



(8.11), respectively. Plotted is the energy multiplied by QyjS^ as a function of Q. The square 



markers show the exact results, which lie very close to the physical branch in each case. 

modes in general. So this is a formidable task. Higher-dimensional analogues, polytopes, 
are also currently the subject of ongoing work. 
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Appendix A: Poisson Resummation Formulas 



We consider the Poisson resummation of the traced cylinder kernel of an arbitrary real 
quadratic form, 

oo 



m,i,...,m„=— oo 



Taking the Fourier transform of the summand of S and using Eq. (1.2) gives 



^ /"OO ^ 

mi,...,m„=— oo 3=^ 



We shift the variables 



Uj — 7- Uj — aj , 



(A2) 



(A3) 
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and diagonalize A 

B.j^UikAkmUlj. (A4) 
A redefinition of tlie integration variables follows, 

Vj = UjkUk , (A5) 

as well as the summation variables, 

qj^UjkrUk. (A6) 

As a result of these transformations, we recognize that the Jacobian of the transformation 
matrix is unity, 

n n 

Hduj^Hdvj. (A7) 

i=i 3=1 

We are now ready to change to hyperspherical coordinates. First, we define 

Rj = VWj ^3 (A8) 

and 

% = (A9) 

which allows us to write 

VjQj^kRcosO. (AlO) 
Effectuating the change of variables gives us 

n n—3 

Y[dvj = |det {B)\-^'^ BJ'-^dRd(t) {smOf-^dO \[ {smO^y ddj . (All) 
The 0-integral produces a 27r and the integrals for the first (n — 3) Oj angles give 

n(/ "'«') = fo^r^rp)- 







3=1 

We are now able to focus on the remaining ^-integral. 



(sin ef-'' 6^'^^*^^'=°^^ de = 7r(^-")/'r((n - 1) /2) {kRf-''^^^ Jin-2)/2 i^nkR) . (A13) 

Jo 

The last integral, the i?-integral, is evaluated rather straightforwardly, 

Jo ''''' - (,2+4,.,^.,^,)(n-.l)/2 ' 
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and putting everything together we obtain: 

2"7r("-l)/2r((n+ l)/2) ^ ^g-2^im,a,^ 



From this result we obtain the following resummed expressions we use in the paper: 



V drj ^ y/a^T* 2n3y/^ ^ ' 

^ ' p,q,r=—oo \ I I \ I I 



oo 

X 

p,q,r=—oo 



E' 



{pya + qy/3 + rV7) 



' p=—oo * p=— oo 

Here the prime means that all positive and negative integers are included in the sum, but 
not the case where all the integers are zero. 



Appendix B: Epstein Zeta Functions 

We define the following Epstein zeta functions: 

oo 

Z3(s;a,6,c)= ^ {a + bm^ + cn'^y , (Bl) 

A;,m,n=— oo 

oo 

Z3h{s;a,b,c) = (-l)'=+"+"(aP + 6m2 + cn2)-^ (B2) 

fc,m,n=— oo 
oo 

Z3c{s]a,b,c) = {-l)''+'^{ae + bm^ + cnY% (B3) 

A;,m,n=— oo 
oo 

Z2b{s;a,b)= (-ir+"(a^' + &^')"'- (B4) 

m,n=— oo 

Here, sums are over all integers, positive, negative, and zero, excluding the single point 
where all are zero. They are summed numerically using Ewald's method [191 |25l |26]. A few 
specific values needed for calculations: 

^3(2; 1, 1, 1) = 16.5323159598, (B5) 
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^3fc(2; 1,1,1 



3.8631638072, 



(B6) 



^3c(2; 1,1,1 



1.8973804658, 



(B7) 



^2^(3/2; 1,2) 



1.9367356117, 



(B8) 



^26(3/2; 1,3) 



1.8390292892. 



(B9) 



The Dirichlet L-series are defined as Lk{s) = Yl^=iXk{n)n where Xfc is the number- 
theoretic character [19]. The Dirichlet beta function, also known as L_4, is usually defined 

as/3(5) = Er=o(-ir(2n + l)~'^- 
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